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We study the director field around a spherical particle immersed in a uniformly aligned nematic 
liquid crystal and assume that the molecules prefer a homeotropic orientation at the surface of the 
particle. Three structures are possible: a dipole, a Saturn-ring, and a surface-ring configuration, 
which we investigate by numerically minimizing the Frank free energy supplemented by a magnetic- 
field and a surface term. In the dipole configuration, which is the absolutely stable structure 
for micron-size particles and sufficiently strong surface anchoring, a twist transition is found and 
analyzed. We show that a transition from the dipole to the Saturn ring configuration is induced 
by either decreasing the particle size or by applying a magnetic field. The effect of metastability 
and the occurence of hysteresis in connection with a magnetic field are discussed. The surface-ring 
configuration appears when the surface-anchoring strength W is reduced. It is also favored by a large 
saddle-splay constant K 24. A comparison with recent experiments p"o| , [i"3[ gives a lower bound for 
W, i.e., W > 0.6 erg/cm 2 for the interface of water and pentylcyanobiphenyl (5CB) in the presence 
of the surfactant sodium dodecyl sulfate. 

PACS numbers: 77.84.Nh, 61.30.Cz, 61.30.Jf 



I. INTRODUCTION 

Dispersions of particles in a host fluid are part of our 
everyday life and an important matter of state. Since 
they appear in food (e.g. milk), drugs, paints, and ink, 
they are of considerable technological importance. In col- 
loidal suspensions the particles are solid, whereas emul- 
sions are dispersions of liquid droplets coated with a 
surfactant. There exists a whole zoo of interactions be- 
tween the particles whose delicate balance determines the 
stability of a dispersion. Besides the conventional van- 
der-Waals, screened Coulombic, and steric interactions 
0, fluctuation- induced Casimir forces |^,|| and deple- 
tion forces have attracted a lot of interest. 

When the particles are dispersed in a nematic liquid 
crystal, i.e., an anisotropic fluid, where elongated organic 
molecules are aligned on average along a common di- 
rection called director, additional long-range forces due 
to elastic deformations of the director field are induced. 
These forces are of dipolar or quadrupolar type depend- 
ing on the symmetry of the director configuration around 
the particle @-[ll]] , and they were confirmed by recent ex- 
periments in inverted nematic emulsions [ p^| , ^2|JT^ ] . On 
the other hand topological point defects in the orienta- 
tional order of the molecules give rise to a short-range 
repulsion 10 13]. Related early work deals with forces 



between walls mediated by an enclosed nematic liquid 
crystal JlJ,|l^|. Surface- induced nematic order in the 
isotropic phase close to the isotropic-nematic phase tran- 
sition leads to a short-range interaction as reported in 
most recent experimental and theoretical work [|l6] p"9[ . 
Even Casimir forces arising from fluctuations in the liquid 
crystalline order parameter have been calculated [g0| |2[ . 

Before treating the interaction between particles it is 
necessary to understand how one single particle behaves 
in a nematic environment. This article deals with the 



possible director configurations around a spherical par- 
ticle, which is placed into a uniformly aligned nematic 
liquid crystal, and which prefers a homeotropic, i.e., per- 
pendicular anchoring of the molecules at its surface. This 
is already a complicated problem. Since it cannot be 
solved analytically without employing an ansatz func- 
tion for the director configuration and further approx- 
imations, we will handle it by numerical methods. 

Poulin et al. showed experimentally that in inverted ne- 
matic emulsions, where surfactant-coated water droplets 
are dispersed in a nematic liquid crystal, a director field 
configuration of dipolar symmetry occurs (see Fig. |l|) 
p0| , p3[ . The water droplet and its companion hyperbolic 
point defect form a tightly bound object which we call 
dipole for short. A similar observation at a nematic- 
isotropic interface was made by R. B. Meyer much ear- 
lier (23|. Both the droplet and the defect carry a topo- 
logical charge +1, which "add up" to the total charge 
of the dipole |24||. Theoretically it was described with 
the help of ansatz functions that were motivated by an 
electrostatic analog ]lo|Jri| . Terentjev et al. introduced 
the Saturn-ring configuration with quadrupolar symme- 
try where a —1/2 disclination ring surrounds the sphere 
at the equator (see Fig. [I]). It was investigated by both 
analytic and numerical methods p5| , |26| . By shrinking 
the disclination ring to the topologically equivalent hy- 
perbolic point defect, the Saturn ring can be continuously 
transformed into the dipole configuration jnj . With the 
help of an ansatz function that describes such a transfor- 
mation it was conjectured that for sufncently small parti- 
cles the Saturn ring should be more stable than the dipole 
jy}. For a finite anchoring strength of the molecules at 
the surface a third structure occurs, which we also illus- 
trate in Fig. [I]. We call it the surface-ring configuration. 
Depending on the anchoring strength W there exists ei- 
ther a —1/2 disclination ring sitting directly at the sur- 
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face, or, for smaller W p6| , the director field is smooth 
everywhere, and a ring of tangentially oriented molecules 
is located at the equator of the sphere. By means of a 
Monte-Carlo simulation Ruhwandl and Terentjev showed 
that for sufficently small anchoring strength W the sur- 
face ring is the preferred configuration |2j] . 

In this article we give a full account of the three direc- 
tor configurations in Fig. |l| by numerically minimizing the 
Frank free energy supplemented by a magnetic-field and 
a surface term ^8|]. We go beyond the one-constant ap- 
proximation, generally used in the work cited above, and 
include the saddle-splay term of the Frank free energy. 
Furthermore, the results of the analytic approach based 
on the use of ansatz functions are checked. In particular, 
we investigate the dipole configuration, which undergoes 
a twist transition. Then, we show in detail that the tran- 
sition from the dipole to the Saturn ring configuration 
can be either achieved by decreasing the particle size or 
by applying a magnetic field. The role of metastabil- 
ity is discussed. Finally, the surface ring is considered, 
and the special role of the saddle-splay free energy for its 
occurence is pointed out. Lower bounds for the surface- 
anchoring strength W are given. All these results are 
presented in Sec. III. In Sec. |n] we define the geometry 
of our problem, write down the reduced free energy, and 
explain the numerical method to minimize it. 
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FIG. 1. A spherical particle with a preferred homeotropic 
anchoring at its surface that is placed into a uniformly aligned 
nematic liquid crystal exhibits three possible structures: the 
dipole configuration where the particle is accompanied by a 
hyperbolic hedgehog, the Saturn-ring configuration where the 
particle is surrounded by a —1/2 disclination ring at the equa- 
tor, and the surface-ring configuration. 




FIG. 2. Illustration of the modified spherical coordinates 
(p = l/r 2 ,8). At each point the local coordinate basis 
(e r , ee , e<p) of the standard spherical coordinate system is at- 
tached. The integration area is given by < p < 1. 



II. GEOMETRY, FREE ENERGY, AND 
NUMERICAL METHOD 

The director field around a spherical particle follows 
from the minimization of the Frank free energy supple- 
mented by a magnetic-field and a surface term. In this 
section we describe our coordinate system, review the 
free energy, and give some numerical details. 



A. Geometry 

The region outside the spherical particle with radius a 
is infinitely extended. We use a modified spherical coor- 
dinate system with a radial coordinate p — 1/r 2 where 
r is the distance of a space point from the center of the 
particle measured in units of a. The exponent 2 is mo- 
tivated by the far- field of the dipol configuration JlO|,|ll| . 
Such a transformation has two advantages. The exte- 
rior of the particle is mapped into a finite region, i.e., 
the interior of the unit sphere (p < 1). Furthermore, 
equally spaced grid points along the coordinate p result 
in a mesh size in real space which is small close to the 
surface of the particle. In this area the director field is 
strongly varying, and hence a good resolution for the nu- 
merical calculation is needed. On the other hand, the 
mesh size is large far away from the sphere where the 
director field is nearly homogeneous. Since our system 
is axially symmetric, the director field only depends on 
p and the polar angle 9 as illustrated in Fig. ||. The 
symmetry axis corresponds to the z axis. At each point 
(p, 9) we attach the local coordinate basis (e r ,e0,e0) of 
the standard spherical coordinate system and express the 
director in this basis: 



n(p, 9) = n r (p, 9)e r + n e (p, 9)e e + n^p, 9)e 4 



(1) 
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Since the director is a unit vector, we write 

rt. r = cos0 , ng = sin cos $ , n^ — sin sin $ , (2) 

where 0(p, 0) and <&(/?, 0) denote, respectively, the tilt 
and the twist angle. At the surface of the particle we 
allow the director to rotate away from the preferred ra- 
dial direction by introducing a surface free energy (see 
next subsection). At infinity, i.e., at p = n always 
points along the z axis. For completeness we note that 
the differentials dr and dp are connected via 

dr = -dp/(2p 3/2 ) . (3) 



B. Free energy 

The free energy which we will minimize consists of bulk 
and surface terms: 



F n = J d 3 x(f el + f 24 + fn)+J dSfs 



(4) 



with the free energy densities: 
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fel = 2^i(V • n ) + ^(n • V x n) 
+K 3 [n x (V x n)] 2 } 



(5a) 



fu = —zK^V ■ [n(V • n) + n x (V x n)] (5b) 



.fe = --A X [(n.H) 2 -i/ 2 ] 



fs = -W[l - (n • Of] 



(5c) 



(5d) 



that favors an alignment of the director n parallel to 
H. Since we calculate the magnetic free energy of the 
infinitely extended region around the sphere, we use the 
magnetic free energy of a completely aligned director field 
as a reference point to avoid infinities. As a result the 



term —AxH /2 in Eq. 



occurs. Finally, we employ 



the surface free energy of Rapini-Papoular to take into 
account the interaction of the director with the bound- 
aries. In Eq. (5d) the unit vector v denotes some pre- 
ferred orientation of the director at the surface, and W 
is the coupling constant. It varies in the range 1CP 4 — 1 
erg/cm 2 as reviewed by Blinov et al. |3(]]. However, the 
authors do not specify W for the interface of water and 
the liquid-crystalline phase of 5CB in the presence of the 
surfactant sodium dodecyl sulfate, which was used i n the 
experiment by Poulin et al. IToJT^]. In subsection III C 



we will give a lower bound for W for such an interface. 

For the numerical minimization a reduced form of the 
free energy of Eq. (||) is used. We introduce the energy 
unit nK 3 a and refer all lengths to the radius a of the 
spherical particle. Furthermore, we employ the modified 
spherical coordinates (/?,#), take into account the axial 
symmetry of our system, and, finally, arrive at the re- 
duced free energy 



F n = F n /{7rK 3 a) 



d6dp S -^(f el + f 2i + f H ) 



(7) 

d9 sin Of s . 



where 



/ ei = ^i(V-n) 2 +X 2 (n.Vxn) 2 
+ [n x (V x n)] 2 } 



f 24 = -^K 24 V ■ [n(V • n) + n x (V x n)] 



(8b) 



The Oseen-Zocher-Frank free energy density / e ; + f 24 
describes elastic distortions of the director field n(r), 
where Kx, K 2 , K 3 , and K 2 4 denote, respectively, the 
splay, twist, bend, and saddle-splay elastic constants. 
The saddle-splay term is a pure divergence; it, there- 
fore, can be transformed into integrals over all surfaces 
of the system. A Cauchy relation for K 24 follows from 
the Maier-Saupe molecular approach |29|: 



#24 = {Kxx + K 22 )/2 



(6) 



(There is also the possibility of another surface term with 
a free energy density Kx 3 V ■ (nV ■ n), which we will not 
consider in this paper.) Eq. ( |5c| ) couples the director 
to an external magnetic field H, where A% = \\\ ~ X± 
stands for the magnetic anisotropy. The symbols X|| and 
X± denote the magnetic susceptibilities for a magnetic 
field applied, respectively, parallel or perpendicular to 
the director. In this paper we consider a positive A% 



Ih = -2^h 2 [K C0S( 



ng smt 



ft 



(8c) 



(8d) 



In Eqs. ( pa[ ) and ( |8b[ ) the coefficients K\ = Kx/K 3 , K 2 = 
K 2 /K 3 , and K 2 4 — K 2i /K 3 denote, respectively, the 
reduced splay, twist, and saddle-splay elastic constants. 
The single contributions to the Frank free energy density 
in our modified spherical coordinates are rather lengthy, 
and we refer the reader to Appendix ^] for the detailed 
form. We always apply the magnetic field H along the 
symmetry axis of our system, which coincides with the z 
axis (see Fig. ||) . Inserting H = He z with e 2 = cos 9e r — 
sin^ee into Eq. (|5c[) we obtain the magnetic free energy 
density of Eq. (|8c[). The strength of the field is given via 
the reduced magnetic coherence length: 



£ H = ^K 3 /(A x H2)/a 



(9) 
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It indicates the distance in the bulk which is needed to 
orient the director along the applied field when, e.g., 
through boundary conditions a different preferred orien- 
tation of the molecules exists [p8[| . The length tends to 
infinity for H — ► 0. The surface term of Eq. (8d) fol- 
lows from Eq. (5d) by choosing u = e r , i.e., a preferred 



radial anchoring of the molecules at the surface of the 
suspended particle. The reduced extrapolation length 



is = K 3 /{Wa) 



(10) 



signifies the strength of the anchoring. It compares the 
Frank free energy of the bulk, which is proportional to 
Ksa, to the surface energy, which scales with Wa 2 . At 
strong anchoring, i.e., for £5 <C 1 the energy to rotate 
the director away from its preferred direction i> at the 
whole surface would be much larger than the bulk en- 
ergy. Therefore, it is preferable for the system when the 
director points along is. However, n can deviate from i> 
in t he are a £sa 2 . This will explain the result in subsec- 
tion III C where we show that the surface ring already 
appears at strong surface coupling. Rigid anchoring is 
realized for £s — > 0, and (s>l means weak anchoring, 
where the influence of the surface is minor. 

To calculate the free energy F n numerically we trans- 
form the saddle-splay term into a surface integral: 

F 2i = K24 / d8sin6e r ■ [n(V • n) + n x (V x n)] . 

J p=i 

(11) 

The detailed form follows from Appendix [a|. In perform- 
ing Gauss' theorem we only include the surface of the 
sphere. There is also a contribution from the surface of 
the core of a possible disclination ring. However, since 
we can only describe the energy of the core in an approx- 
imate way (see next paragraph), we will skip this term. 
For rigid homeotropic anchoring (n r = 1, rig = = 0) 
the free energy of the saddle-splay term is 

F 2i = 4if 24 . (12) 

To arrive at Eq. ( |l2] ) we used dn r /dp = 0, which is valid 
because of the normalization of the director. 

The free energy of a ±^ disclination ring is taken into 
account by the line energy Fd of a ±i disclination. In 
the one-constant approximation [K\ = K 2 = K$ = K) 
and in reduced units it reads Ik5l| 



1 4 1 2 rv 



(13) 



The first term denotes the line energy of the core with a 
core radius r c which in absolute lengths is of the order 
of 10 nm. The second term stands for the elastic energy 
of the disclination line where R is its radial extent [see 
Fig. |]b)]. In the general case {K\ ^ K 2 7^ K 3 ) an ana- 
lytic expression for the elastic energy does not exist. We 



only use a rough approximation for the core energy F c 
by averaging over the Frank constants: 



F c = F c /K 3 = -(K 1 + K 2 + l)/3 



(14) 



A more quantitative description of the free energy of 
a disclination has to start from the Landau-Ginzburg- 
de Gennes free energy with the full alignment tensor Q 

us 

We finish this subsection with an important remark 
about length scales. In the free energy of Eq. ((?]) all 
lengths are given relative to the particle radius a, which 
appears in the absolute energy unit irK^a only. This 
would suggest that the director configuration does not 
depend on the particle size. However, with the core ra- 
dius r c a second length scale is introduced, which in ab- 
solute lengths is always of the order of 10 nm [|TJ. On 
the other hand, in the next subsection we explain that 
r c in units of a is a parameter of our numerical method 
whose lower bound is determined by the mesh size of the 
grid. In the discussions of Sec. Ill we want to make a con- 
nection to experiments. Therefore, we do not give this 
dimensionless r c , but use an absolute core size of 10 nm 
to calculate the corresponding particle radius a: 

a = 10 nm/r c . 



(15) 



In subsection III B 1 we will study how the configuration 
around a spherical particle depends on a. 



C. Minimization and numerical details 

The fields of the tilt [Q(p, 9)] and the twist [$(p, 9)} an- 
gle follow from a minimization of the reduced free energy 
F n [0,<i>] of Eq. (Q). The corresponding Euler-Lagrange 
equations are equivalent to the functional derivatives of 
F a [G,$]: 



SF n _ SFndjH 

(59 
6F n 



SF n drii 



(5$ 



<9$ 



= 







(16) 
(17) 



where i stands for r, 9, and <fi and Einstein's summation 
convention over repeated indices is used. We have em- 
ployed a chain rule to arrive at the Euler-Lagrange equa- 
tions for <d(p, 9) and Q(p, 9), which is generally valid for 
functional derivatives as shown in Appendix [b]. Perform- 
ing the variation of the free energy F n we arrive at the 
Euler-Lagrange equations in the bulk, 
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d \ 
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dn it3 j 
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(fel + f H ) 
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(18) 
(19) 
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FIG. 3. a) Coordinates (rd,6d) for a -1/2 disclination ring 
with a general position around the spherical particle; b) Cross 
section of the —1/2 disclination ring with a radial extension 
R. 



and at the surface, 



dfs 

dm 

dfs 
dm 



= 



8 



On 



hP 



dm 
d<S> 



(20) 
(21) 



where m ,j stands for dm/dxj . The Euler-Lagrange equa- 
tions are calculated with the help of the algebraic pro- 
gram Maple and are then imported into a Fortran pro- 
gram where they are usually solved on a rectangular 
33 x 129 grid in the (p, 0) space. In Eq. (33) of ref. O] 
an analytical form of the director configuration of a 
disclination ring around a spherical particle is given. The 
position of this ring is determined by the radial (r d ) and 
the angular (6 d ) coordinates [see Fig. ||a)]. We take this 
director configuration and let it relax via the standard 
Newton-Gauss-Seidel method p3j. In the dipole config- 
uration (6d = 0, 7r) the hyperbolic point defect moves 
along the z axis in its local minimum when the numeri- 
cal relaxation is performed. However, a disclination ring 
(6 d 7^ 0, 7r) basically stays at the positon where we place 
it. We use this fact to investigate the free energy as a 
function of r d and 9 d which gives an instructive insight 
into potential barriers for a transition between the dipole 
and the Saturn-ring configuration. 

The free energy F n of the director field follows from 
a numerical integration. This procedure assigns some 
energy to the disclination ring which certainly is not cor- 
rect. To obtain a more accurate value for the total free 
energy F we use the formula 



F = F n — F n \ + F c/d x 2nr d sin9 d 



(22) 



The quantity -F n | torus denotes the numerically calcu- 
lated free energy of a toroidal region of cross section ttR 2 
around the disclination ring [see Fig. [|b)] . Its volume is 
ttR 2 x 2Trr d sm0 d , where the coordinates (r d ,8 d ) of the 
ring are determined by searching for the maximum of 
the local free energy density f el . The value -F n | t is 
replaced by the last term on the right-hand side of Eq. 



( ^2| ) , which provides the correct free energy according to 
Eqs. ( |l3| ) or (|l4|). To find out how large the cross sec- 
tion ttR 2 of the cut torus has to be, we employed the last 
formula and the line energy of Eq. ( |l3|) for constant r c 
and varying R. To be consistent, F should not depend 
on R. Within an error of less than 1 % this is the case 
if ttR 2 is equal or larger than 3ApA#/2 where Ap and 
AO are the lattice constants of our grid. To study the 
transition between the dipole and the Saturn ring as a 
function of the particle size we choose ttR 2 = 25ApA9/2, 
employ Eq. ( |l3| ) for different r c , and calculate a from Eq. 
(|l5). In all other cases we set R = r c , determine r c from 
7rrj = 3ApA#/2, and use the core energy of Eq. (jlj). 

The radial extension of the core of a point defect is also 
of the order of 10 nm [^5| , and its free energy is approxi- 
mated by A3 x 10 nm. As we show in the following section 
the free energy of the dipole amounts to around 107rA3a. 
Since we consider particle radii larger than 100 nm, the 
contribution from the energy of the point defect is al- 
ways smaller than 1 %. This is beyond our numerical 
accuracy, and therefore no energetical correction for the 
point defect was included. 

The discussion in the following section always uses the 
one-constant approximation or the nematic liquid crys- 
tal pentylcyanobiphenyl (5CB) with the bend elastic con- 
stant A3 = 0.53 x 10~ 6 dyn and the reduced splay and 
twist elastic constants Ai = 0.79 and A 2 = 0.43, respec- 
tively. 



III. RESULTS AND DISCUSSION 

In this section we present the results of our numerical 
investigation. We first address a twist transition in the 
dipole configuration. Then we discuss a transition from 
the dipole to the Saturn ring which is induced either by 
decreasing the particle size or by applying a magnetic 
field. Finally we illustrate that the surface-ring config- 
uration appears when the surface-anchoring strength is 
lowered. 



A. Twist transition of the dipole configuration 

In Fig. H we plot the director field of the dipole config- 
uration for the one-constant approximation. A magnetic 
field is not applied and the directors are rigidly anchored 
at the surface. The dot indicates the location of the hy- 
perbolic hedgehog. For its distance r d from the center 
of the sphere we find r d — 1.26 ± 0.02, where the mesh 
size of the grid determines the uncertainty in r d . Our 
result is in excellent agreement with Ref. Hj. In this 
article the dipole was described via an ansatz function. 
However, Ruhwandl and Terentjev using a Monte-Carlo 
minimization report a somewhat smaller value for r d [pj. 

Fig. |^ presents the distance r d as a function of the re- 
duced splay (Ai) and twist (Ai) constants. In front 
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FIG. 4. The numerically determined director field of the 
dipole configuration for the one-constant approximation, zero 
magnetic field, and rigid surface anchoring. The location of 
the hyperbolic hedgehog is indicated by a dot. The directors 
lie in the drawing plane. 

of the thick line is basically constant. Beyond the line 
rd starts to grow which indicates a structural change in 
the director field illustrated in the nail picture of Fig. [fj. 
Around the hyperbolic hedgehog the directors develop 
a non-zero azimuthal component introducing a twist 
into the dipole. It should be visible under a polarizing 
microscope when the dipole is viewed along its symmetry 
axis. 

In Fig. |?] we draw a phase diagram for the twist transi- 
tion. As expected, it occurs when K\ increases or when 
K2 decreases, i.e., when a twist deformation costs less 
energy than a splay distortion. The open circles are nu- 




FIG. 5. The distance rd of the hyperbolic hedgehog from 
the center of the sphere as a function of the reduced splay 
(Ki) and twist (K2) constants. 
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FIG. 6. Nail picture of a close-up of the twisted dipole 
configuration. Around the hyperbolic hedgehog the directors 
are tilted out of the drawing plane. The length of the nail is 
proportional to the projection of the director on the drawing 
plane. The head of the nail is below the plane. 

merical results for the transition line which can be well 
fitted by the straight line K2 ~ K\ — 0.04. Interest- 
ingly, the small offset 0.04 means that K 3 does not play 
an important role. Typical calamatic liquid crystals like 
MBBA, 5CB, and PAA should show the twisted dipole 
configuration. 

Since the twist transition breaks the mirror symmetry 
of the dipole, which then becomes a chiral object, we 
describe it by a Landau expansion of the free energy: 

F = F + a(K x ,K 2 )[n^] 2 + c[n^] 4 . (23) 

With the maximum azimuthal component n™ ax we have 
introduced a simple order parameter. For symmetry rea- 
sons only even powers of n™ ax are allowed. The phase 

transition line is determined by 0(^1,^2) = 0. Accord- 
ing to Eq. ( p3| ) we expect a power-law dependence of the 
order parameter with the exponent 1 /2 in the twist region 




*1 

FIG. 7. Phase diagram of the twist transition as a function 
of the reduced splay (Ki) and twist (K2) constants. A full 
explanation is given in the text. 
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FIG. 8. Log-log plot of the order parameter n™ ax versus 
Ki close to the twist transition for Ki = 0.4; o . . . numerical 
values, — ... fit by a straight line. 



close to the phase transition. To test this idea we choose 
a constant K% and determine n™ ax for varying K\. As 
the log-log plot in Fig. || illustrates, when approching the 
phase transition, the order parameter obeys the expected 
power law: 



(K 1 - 0.4372) 



1/2 



K 2 = 0.4 



(24) 



B. Dipole versus Saturn ring 

There are two possibilities to induce a transition from 
the dipole to the Saturn-ring configuration; either by re- 
ducing the particle size or by applying, e.g., a magnetic 
field. We always assume rigid anchoring in this subsec- 
tion, set K24 — 0, and start with the first point. 



1. Effect of particle size 

In Fig. [| we plot the reduced free energy F as a func- 
tion of the angular coordinate 9d of the disclination ring. 
For constant 9d the free energy F was chosen as the min- 
imum over the radial coordinate r<j. The particle radius 
a is the parameter of the curves and the one-constant 
approximation is employed. Recall that 9d = n/2 and 
9d = 7r correspond, respectively to the Saturn-ring or 
the dipole configuration. Clearly, for small particle sizes 
(a = 180 nm) the Saturn ring is the absolutely stable 
configuration and the dipole enjoys some metastability. 
However, thermal fluctuations cannot induce a transition 
to the dipole since the potential barriers are much higher 
than the thermal energy k^T . E.g., a barrier of OAirK^a 
corresponds to 1000 k B T (T = 300 K, a = 1/jm). At 
a ss 270 nm the dipole assumes the global minimum of 
the free energy, and finally the Saturn ring becomes abso- 
lutely unstable at a w 720 nm. The scenario agrees with 
the findings of Ref. Ml. Furthermore, we stress that 
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FIG. 9. The reduced free energy F as a function of the 
angular coordinate 6d- The parameter of the curves is the 
particle size a. Further parameters are indicated in the inset. 




the particle sizes were calculated from Eq. (15) with the 



choice of 10 nm as the real core size, and that our results 
depend on the line energy ( ji~3| ) of the disclination. 

The preferred radial coordinate of the disclination 
ring as a function of 9d is presented in Fig. [l(]. As long 
as the ring is open Td does not depend on 9d within an 
error of ±0.01. Only in the region where it closes down to 
the hyperbolic hedgehog does Td increase sharply. The 
figure also illustrates that the ring sits closer to larger 
particles. The radial position of = 1.10 for 720nm 
particles agrees very well with Refs. jn| and ||. 



2. Effect of a magnetic field 

A magnetic field applied along the symmetry axis of 
the dipole can induce a transition to the Saturn-ring 
configuration. This can be understood from a simple 
back-of-the-envelope calculation. Let us consider high 
magnetic fields, i.e., magnetic coherence lengths much 
smaller than the particle size a, which in our reduced 
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FIG. 10. The preferred radial coordinate of the discli- 
nation ring in units of a as a function of 0d for two particle 
sizes. Further parameters are indicated in the inset. 
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FIG. 11. The reduced free energy F as a function of the 
angular coordinate 6d- The parameter of the curves is the 
inverse reduced magnetic coherence length fTj 1 . Further pa- 
rameters are indicated in the inset. 

units means: <C 1. The directors are basically aligned 
along the magnetic field. In the dipole configuration the 
director field close to the hyperbolic hedgehog cannot 
change its topology. The field lines are "compressed" 
along the z direction, and high densities of the elastic and 
magnetic free energies occur in a region of thickness 
Since the field lines have to bend around the sphere the 
cross section of the region, in units of a 2 , is of the order 
of 1, and its volume is proportional to The Frank 
free energy density is of the order of l/£fj, and therefore 
the elastic free energy, in reduced units, scales with 1/£h- 
The same holds for the magnetic free energy. In the case 
of the Saturn- ring configuration high free energy densities 
occur in a toroidal region of cross section oc £L around 
the disclination ring. Hence, the volume scales with £jj 
and the total free energy is of the order of 1, i.e., a factor 
1/£h smaller than for the dipole. 

Fig. [n] presents a calculation for a particle size of 
a = 0.5 /xm and the liquid crystal 5CB. We plot the re- 
duced free energy as a function of 9d for different mag- 
netic field strengths given in units of the inverse reduced 
coherence length ^ . Without a field (£7^ = 0) the 
dipole is the energetically preferred configuration. The 
Saturn ring shows metastability. A thermally induced 
transition between both states cannot happen because of 
the high potential barrier. At a field strength Q 1 = 0.33 
the Saturn ring becomes the stable configuration. How- 
ever, there will be no transition until the dipole looses its 
metastability at a field strength t^ 1 — 3.3, which is only 
indicated by an arrow in Fig. [ll]. Once the system has 
changed to the Saturn ring it will stay there even for zero 
magnetic field. Fig. |l2]a) schematically illustrates how a 
dipole can be transformed into a Saturn ring with the 
help of a magnetic field. If the Saturn ring is unstable 
at zero field, a hysteresis occurs [see Fig. |l2|b)]. Com- 
ing from high magnetic fields the Saturn ring looses its 
metastability at H t \ and a transition back to the dipole 
takes place. In Fig. ^ we showed that the second situa- 
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FIG. 12. a) The Saturn ring is metastable at H = 0. The 
dipole can be transformed into the Saturn ring by increas- 
ing the magnetic field H beyond Ht2 where the dipole looses 
its metastability. Turning off the field the Saturn ring stays, 
b) The Saturn ring is unstable at H = 0. When the mag- 
netic field is decreased from values above Ht2, the Saturn 
ring shrinks back to the dipole at Hti where the Saturn ring 
looses its metastability. A hysteresis occurs. 

tion is realized for particles larger than 720 nm. We also 
performed calculations for a particle size of 1/xm and the 
liquid crystal 5CB and still find the Saturn ring to be 
metastable at zero field in contrast to the result of the 
one-constant approximation. 

Finally, in Fig. |l3| we plot the reduced free energy ver- 
sus the applied magnetic field for different particle sizes. 
The energy of the dipole (dashed line) does not depend 
on a. However, the field strength where the dipole be- 




FIG. 13. The reduced free energy F as a function of the 
inverse coherence length fjj for the dipole and the Saturn 
ring. The parameter of the curves is the particle size a. Fur- 
ther parameters are indicated in the inset. 
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comes unstable should depend on the particle size. The 
dot indicates this strength for a — 0.5 fim. The energy 
curves of the Saturn ring are merely shifted by a con- 
stant amount when the particle size a is changed. For 
a = 0.25 /xm the Saturn ring is the stable configuration. 
At the intersection points of full and dashed curves the 
absolute stability changes from the dipole to the Saturn 
ring. These points seem to occur at higher fields when 
the particle size is increased. However, if we take into ac- 
count that ^Jj 1 oc H is given relative to there is not 
much variation in the absolute field strength H between 
the two particle sizes of 0.5 and 1 /im. Larger particles 
could not be investigated because they would have re- 
quired smaller mesh sizes of the grid. 



C. Influence of finite surface anchoring 



viates from the homeotropic alignment (see inset of Fig. 
|l5| ) then a simple energetical estimate allows AO to be 
of the order of £s. It is interesting that the transition 
point shifts to higher anchoring strengths, i.e, decreasing 
£s when K24 is increased. Obviously the saddle-splay 
term favors the surface-ring configuration. To check this 
conclusion we plot in Fig. |l5| the saddle-splay free energy 
F24 versus £g. The horizontal lines belong to the dipole. 
They correspond to the energy which one expects for a 
rigid homeotropic anchoring at the surface of the sphere 
[see Eq. |l2|)]. In contrast, for the surface-ring configura- 
tion the saddle-splay energy drops sharply. The surface 
ring around the equator of the sphere introduces a "sad- 
dle" in the director field as illustrated in the inset of 
Fig. 15. Such structures are known to be favored by the 



saddle-splay term. We modeled the surface ring with an 
angular width AO by the following director components: 



In the last subsection we investigate the effect of fi- 
nite anchoring on the director field around the spherical 
particle. The saddle-splay term with its reduced constant 
K24 is important now. We always choose a zero magnetic 
field. In Fig. ^ we employ the one-constant approxi- 
mation and plot the free energy versus the surface ex- 
trapolation length £5 for different saddle-splay constants 
K24- Recall that £5 is inversely proportional to the sur- 
face constant W [see Eq. ([To|)]. The straight lines belong 
to the dipole. Then, for decreasing surface anchoring, 
there is a first-order transition to the surface-ring struc- 
ture. We never find the Saturn-ring to be the stable 
configuration although it enjoys some metastability. For 
K24 = the transition takes place at £s « 0.085. This 
value is somewhat smaller than the result obtained by 
Ruhwandl and Terentjev ||. One could wonder that the 
surface ring already occurs at such a strong anchoring like 
£s ~ 0.085 where any deviation from the homeotropic 
anchoring costs a lot of energy. However, if AO is the an- 
gular width of the surface ring where the director de- 



tanh 



tt/2 



A6 



cosh 



-tt/2 
A9 



(25a) 



(25b) 



where Ad <§C tt/2 to ensure that n r = 1 at — 0, n, and 
calculated the saddle-splay energy versus A9 by numer- 
ical integration. The result for K24 = 1 is shown in the 
inset of Fig. [l5]. It fits very well to the full numerical 
calculations and confirms again that a narrow "saddle" 
around the equator can considerably reduce the saddle- 
splay energy. 

For the liquid crystal 5CB we determined the stable 
configuration as a function of K24 and £g. The phase 
diagram is presented in Fig. [l(| With its help we can 
derive a lower bound for the surface constant W at the 
interface of water and 5CB when the surfactant sodium 
dodecyl sulfate is involved. As the experiments by Poulin 




0.05 0.1 0.15 0.2 0.25 0.3 



FIG. 14. The minimum free energy F as a function of the 
surface extrapolation length £s for different K24- A first-order 
phase transition from the dipole to the surface ring occurs. 
Further parameters are indicated in the inset. 
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FIG. 15. The saddle-splay free energy F24 as a function of 
£s for the same curves as in Fig. [u| Inset: F24 versus the 
angular width of the surface ring calculated from the ansatz 
functions in Eqs. ( |25| ) for K24 = 1. 
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et al. clearly demonstrate water droplets dispersed in 
5CB do assume the dipole configuration. From the phase 
diagram we conclude £s < 0.09 as a necessary condition 
for the existence of the dipole. With awl |xm, K3 = 
0.53 x 10~ 6 dyn, and the definition (10) for £g we arrive 
at 



W > 0.6 erg/cm 2 



(26) 



If we assume the validity of the Cauchy-Relation (||) 
which for 5CB gives K24 = 0.61, we conclude that 
W > 1.5 erg/cm 2 . Both values indicate one of the highest 
anchoring strengths ever observed. 




0.1 



0.04 0.06 0.08 

FIG. 16. Phase diagram of the dipole-surface ring transi- 
tion as a function of £s and K24- Further parameters are 
indicated in the inset. 



IV. CONCLUSIONS 



position 8 ( i or is even shrunk to a point at 9^ — 0, ir. The 
surface ring is clearly favored by a large saddle-splay con- 
stant -^24- 

So far convincing experiments on dispersions of spher- 
ical particles in a nematic liquid crystal only exist in the 
case of inverted nematic emulsions [[H)]|iJ,[f| . We hope 
that the summary of our results stimulates further exper- 
iments which try different liquid crystals as a host fluid, 
manipulate the anchoring strength, investigate the effect 
of external fields, and attempt to disperse silica or latex 
spheres p6|-p8|. 
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APPENDIX A: 

For completeness we give the explicit formulas of V • n, 
n • V x n, and n x (V x n) for the director field in Eq. 
(|l|) and the modified spherical coordinates (p, 9) : 

V • n = v / p(2n r — 2pn r:P + cot 9ng + ng t g) (Ala) 



n • V x n = ^fp(2pn s n <i> , p - 2pn < j > ng tP - n^n rfi 

+ cot 6n r ntj) + n r n^ t e) (Alb) 



The purpose of the article was to give a detailed study 
of the director field around a spherical particle and to 
illustrate how it can be manipulated by an external field. 
We clearly find that for large particles and sufficiently 
strong surface anchoring the dipole is the preferred con- 
figuration. For conventional calamitic liquid crystals 
where K% < K\ the dipole should always exhibit a twist 
around the hyperbolic hedgehog. It should not occur 
in discotic liquid crystals where K2 > K\. According to 
our calculations the bend constant K3 only plays a minor 
role for the twist transition. The Saturn ring appears for 
sufficiently small particles. However, the dipole can be 
transformed into the Saturn ring by means of a magnetic 
field if the Saturn ring is metastable at H = 0. Other- 
wise a hysteresis is visible. For the liquid-crystal 5CB 
we find the Saturn ring to be metastable for a particle 
size a = 1 pxa. Increasing the radius a this metastability 
will vanish in analogy with our calculations within the 
one-constant approximation (see Fig. ^). Decreasing the 
surface-anchoring strength W the surface ring configu- 
ration with a quadrupolar symmetry becomes absolutely 
stable. We never find a stable structure with dipolar 
symmetry where the surface ring has a general angular 



n x (V x n) = *Jp[ne(ne - 2png, p - n r ,e) 
+12^(71^ - 2pn0 iP )]e r 
+^/p[n^(cot On^ + n^g) 
-n r (n e - 2png tP - n r ,g)]e g 
-y/p[n r (n<f, - 2pn0 :P ) 
+ng(cot 0710 + 7i0,e)]e^ . (Ale) 

Partial derivatives are indicated by , i (i — p, 9) and 

I - -^ m h < A2 » 

was used. 



APPENDIX B: CHAIN RULE FOR 
FUNCTIONAL DERIVATIVES 

Suppose F[$i] is a functional depending on the fields 
3>i(x) in real space. The functional derivative 8F/6$i(x.) 
is introduced via the Taylor expansion 



10 



(Bl) 

8F indicates the change of the functional if at every space 
point the field $i(x) is changed by the small amount 
<5$i(x). The special choice <5$; = 5ijc5(x — xo) leads 
directly to the definition 



5F 



:= lim 

e->0 



F[$i + %e5(x - x )] - F[$J 



(B2) 



If the field $j depends on other fields i.e., $j = 
$i(0fe), then with 



F[$i(0 fe + 8 jk e8{i[. - x ))] « F[$,(0 fc ) 



+5 jk e5(x-x ) — -] (B3) 



90, 



one obtains immediately the chain rule 
(5F _ ^ <5F <9$; 



(B4) 



from Eq. (B2) 
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